Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 17 April 2013 (MN I^TeX style file v2.2) 

A scaling-relation for disc galaxies: circular- velocity 
gradient versus central surface brightness 



Federico Lelli^*, Filippo Fraternali^'^, Marc Verheijen^ 

^ Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV, Groningen, The Netherlands 
^Department of Physics and Astronomy, University of Bologna, via Berti Pichat 6/2, 4-0127, Bologna, Italy 



m 

o 

(N 

< 

in 



O 
U 

O 



> 

o 
in 

o 
m 



% 



ABSTRACT 

For disc galaxies, a close relation exists between the distribution of light and the 
shape of the rotation curve. We quantify this relation by measuring the inner circular- 
velocity gradient df(y{Q) for spiral and irregular galaxies with high-quality rotation 
curves. We find that diiV{Q) correlates with the central surface brightness /xq over 
more than two orders of magnitude in d[iV{Q) and four orders of magnitudes in jxq. 
This is a scaling-relation for disc galaxies. It shows that the central stellar density of a 
galaxy closely relates to the inner shape of the potential well, also for low-luminosity 
and low-surface-brightness galaxies that are expected to be dominated by dark matter. 
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1 INTRODUCTION 

Scaling-relations are an ideal tool to investigate the struc- 
ture, the formation, and the evolution of galaxies. For disc 
galaxies, the TuUy-Fisher (TF) relation (TuUy & Fisher 
1977) is one of the best-studied scaling laws. It was originally 
proposed as a correlation between the absolute magnitude 
of a galaxy and the width of its global H I line profile. It is 
now clear that the fundamental relation is between the total 
baryonic mass of the galaxy and the circular velocity along 
the flat part of the outer rotation curve (Vflat), thought to 
be set by the dark matter (DM) halo (e.g. McGaugh et al. 
2000; Verheijen 2001; Noordermeer & Verheijen 2007). 

While Vfizx is related to the total dynamical mass of a 
galaxy, the inner shape of the rotation curve provides in- 
formation on the steepness of the potential well. For disc 
galaxies (Sb and later types), the rotation curve is generally 
described by an inner rising part (nearly solid-body) and 
an outer flat part (e.g. Bosma 1981; Begeman 1987; Swa- 
ters et al. 2009). For bulge-dominated galaxies, instead, the 
rotation curve shows a very fast rise in the center, often fol- 
lowed by a decline and the flattening in the outer parts (e.g. 
Casertano & van Gorkom 1991; Noordermeer et al. 2007). 

The relation between the optical properties of a galaxy 
and the shape of its rotation curve has been debated for 
many years (e.g. Rubin et al. 1985; Corradi & Capaccioli 
1990; Persic & Salucci 1991). Several authors pointed out 
that the shape of the luminosity profile and the shape of 
the rotation curve are closely related (Kent 1987; Caser- 
tano & van Gorkom 1991; Broeils 1992; Sancisi 2004; Swa- 
ters et al. 2009). In particular, de Blok & McGaugh (1996) 
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and TuUy & Verheijen (1997) compared the properties of 
two disc galaxies on the same position of the TF relation 
but with different central surface brightness, and found that 
high-surface-brightness (HSB) galaxies have steeply-rising 
rotation curves compared to low-surface-brightness (LSB) 
ones. Thus, for a given total luminosity (or Vflat), an ex- 
ponential light distribution with a shorter scale-length cor- 
responds to a steeper potential well (see also Amorisco & 
Bertin 2010). For HSB spirals, maximum-disc solutions can 
explain the dynamics in the central regions with reasonable 
values of the stellar mass-to- light ratio K'h/ L (e.g. van Al- 
bada & Sancisi 1986; Palunas & Williams 2000), suggesting 
that either baryons dominate the gravitational potential or 
dark matter closely follows the distribution of light. Garrido 
et al. (2005) also found a clear trend between the inner slope 
of the rotation curve and the central surface brightness of 
18 HSB spiral galaxies. For LSB galaxies, maximum-disc so- 
lutions can reproduce the inner parts of the rotation curves, 
but they often require high values of M^/L that cannot be 
explained by stellar population models (e.g. de Blok et al. 
2001; Swaters et al. 2011), leading to the interpretation that 
LSB galaxies are dominated by DM at all radii. Finally, ga- 
laxies with a central "light excess" with respect to the ex- 
ponential disc (e.g. a bulge) show a corresponding "velocity 
excess" in the inner parts of the rotation curve (e.g. Marquez 
& Moles 1999; Swaters et al. 2009). This is usually referred 
to as the "Renzo's rule" (Sancisi 2004): for any feature in 
the luminosity profile of a galaxy there is a corresponding 
feature in the rotation curve, and vice versa. 

In this Letter, we show that the inner circular-velocity 
gradient dnViQ) of a galaxy strongly correlates with the 
central surface brightness /io over more than two orders of 
magnitude in dnViQ) and four orders of magnitude in po, 
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Table 1. Galaxy sample. The entire table is published in the on-line version of the journal. The last column provides references for the 
distance, the surface photometry, and the rotation curves: a) Jacobs et al. (2009); b) Swaters & Balcells (2002); c) Swaters et al. (2009); 
d) Saha et al. (2006); c) Kent (1987); f) Begeman (1987); g) TuUy et al. (2009); h) Kent (1987); i) de Blok et al. (2008); j) TuUy (1988); 
k) Noordermeer & van der Hulst (2007); 1) Noordermeer et al. (2007). See Sect. 2 for details. 
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thereby extending and firmly establishing the correlation 
hinted at by Fig. 8 of Garrido et al. (2005). This is a scahng- 
relation for disc galaxies. We discuss the implications of this 
relation for the stellar and DM properties of galaxies. 



2 DATA ANALYSIS 

2.1 The circular- velocity gradient 

We define the circular-velocity gradient di{l/(0) as the in- 
ner slope of a galaxy rotation curve, i.e. dV/dR for i? — >■ 0. 
dRV{0) can be estimated if the rising part of the rotation 
curve is well-sampled, but this requires high-quality data 
and a careful modelling of the gas kinematics in the inner 
parts. To minimize the uncertainties, we use four samples 
of galaxies with high-quality rotation curves: Swaters et al. 
(2009) (hereafter S09), de Blok et al. (2008) (hereafter dB08 
or THINGS), Verheijen & Sancisi (2001) (hereafter VSOl), 
and Begeman (1987) (hereafter B87). The rotation curves 
were derived using interferometric H I observations and cor- 
rected for beam-smearing effects. We select only galaxies 
viewed at inclination angles i between 40° and 80°, as the 
rotation velocities of face-on discs require a large correction 
for i, while the observed rotation curves of edge-on discs may 
be affected by unseen holes in the central HI distribution. 
We projected each derived rotation curve onto the corre- 
sponding position- velocity diagram to verify that they have 
been properly determined. We exclude the galaxies from S09 
and VSOl with low-quality data {q > 2, see S09). The spiral 
NGC 3521 (from dB08) is close to edge-on in the inner re- 
gions and we neglect the two innermost velocity points. Five 
galaxies are in common between B87 and dB08. We use the 
new rotation curves from THINGS except for two galaxies 
(NGC 2903 and NGC 3198), as the inner parts of their rota- 
tion curves are better traced by B87, who applied a careful 
beam-smearing correction (see Figs. 9 and 12 of dB08). The 
sample of S09 has one object (NGC 2366) in common with 
dB08 and another one (UGC 6446) in common with VSOl; 
we use the rotation curves from S09. In all these cases the 
differences in diiV{0) are, however, within a factor two. The 
final sample comprises 63 galaxies with morphological types 
ranging from Sab to Sd/Im. 

To this high-quality sample, we add 11 rotation curves 
of SO/Sa galaxies (with 40° ^ i ^ 80°) from Noordermeer 
et al. (2007) (hereafter N07). Since early- type galaxies usu- 
ally lack HI emission in their central regions, the rotation 
curves were derived combining Ha long-slit spectroscopy (for 



the inner parts) with H I observations (for the outer parts) . 
We exclude UGC 12043 as the Ha observations of this galaxy 
have very low velocity resolution. The values of dRV{0) for 
early-type galaxies are more uncertain than those for late- 
type galaxies. 

We derive d_RV(0) by fitting the inner rotation curve 
with a polynomial function of the form 



V{R) = Y,an xR" 



(1) 



and consider the linear term ai = lim_R_>o dV/dR = dRy(O). 
The fit is error-weighted and constrained to pass through 
V = Q At R = Q. The value of ai depends on i) the radial 
range used in the fit, and ii) the order of the polynomial m. 
We define _R9o as the radius where the rotation curve reaches 
90% of its maximum velocity, and fit only the points within 
7?90. This choice allows us to maximize the number of points 
along the rising part of the rotation curve without including 
points along the fiat part. Rotation curves with less then 3 
points within Rgo are excluded, as they are not well-resolved 
in the inner parts. 16 galaxies from the high-quality sam- 
ple and 3 galaxies from N07 are excluded by this criterion, 
thus the high-quality and total samples reduce to 47 and 55 
objects, respectively. For a pure exponential disc with scale- 
length 7?d, ^90 — 1.2i?d- Thus, the first fitted point of the ro- 
tation curve is typically at i? < 0.4i?d and the derived value 
of dR'l^(O) is representative of the innermost galaxy regions 
that are accessible by the available rotation curves. To derive 
the best-fit order of the polynomial, we proceed as follows. 
We start with a linear fit (m = 1) and progressively increase 
m until the reduced-x^ ixi) approaches 1. In practice, we 
minimize the function Px(x^i i') — 0.5, where Px{x^\ v) is the 
integral probability of x^ and v is the number of degrees of 
freedom; the procedure is halted in case x^ would drop below 
1. Visual inspection showed that this method works better 
than the F-test (e.g. Bevington & Robinson 2003), that in 
some cases returns high values of m and thus increases the 
number of free parameters in the fit. 

We test our automatic procedure on a set of model ro- 
tation curves, calculated by summing the contributions of 
a disc, a bulge, and a DM halo. We add typical errors to 
the velocity points (~5 km s~^) and try several spatial sam- 
plings. We find that, even if the rotation curve is poorly- 
sampled (~5 points within i?9o), the actual value of da 1/(0) 
can be recovered with a error of ~30%. However, if the ro- 
tation curve has an inner "bump" (due to a compact bulge), 
d-R,V{Q) may be under-estimated by a factor ~2. 

Figure 1 shows the results for four representative ga- 
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Figure 1. Results of the polynomial- fit for four representative ga- 
laxies. Filled circles show the points of the rotation curve within 
Rgo, while open circles show the points excluded in the fit. The 
solid, red line shows the fitted polynomial function, while the 
dashed, blue line shows its linear term. The order m of the poly- 
nomial is indicated. See Sect. 2.1 for details. 



laxies that require polynomial fits of difTerent orders. Late- 
type galaxies (Sb to Im) are generally well-fitted by poly- 
nomials with m, = 1 (e.g. UGC 7559) or m = 2 (e.g. 
NGC 3198), but several cases do require tti ^ 3 (e.g. 
NGC 5055). Early-type galaxies (SO/Sa) often require high- 
order polynomials (m ^ 4, e.g. UGC 11670), as their ro- 
tation curves may have complex shapes characterized by a 
steeply-rising part followed by a decline and a second rise. 
For some bulge-dominated galaxies from N07, the value of 
dRV{Q) is rather uncertain, since there may be no data 
points in the inner radial range where the linear term a\ 
is representative of the rotation curve (see Fig. 1, bottom- 
right). Table 1 provides the fit results for the galaxies in 
Fig. 1. The results for the entire galaxy sample are provided 
in the on-line version of this Table. 

The error Sd^vio) on dijV^(O) is estimated as 



'^-^(o)^) +(^-^(o)|) (2) 



Sd^v{o) = \l5h + 



where 5ax is the nominal error on the fitted linear term ai, 
5i is the error on the disc inclination i, and (5d is the error 
on the galaxy distance D. fo typically gives a negligible 
contribution for galaxies with distances derived using the tip 
of the red giant branch (TRGB) and/or Cepheids (Ceph), 
whereas it can dominate the error budget for galaxies with 
distances estimated from the TF relation or the Hubble fiow. 



2.2 The central surface brightness 

For the high-quality sample of disc-dominated galaxies, we 
consider two ways to estimate the central surface brightness: 
i) the disc central surface brightness /.td, obtained from an 
exponential fit to the outer parts of the luminosity profile; 



and ii) the observed central surface brightness /xq, obtained 
from a linear extrapolation of the luminosity profile in the 
inner few arcseconds to i? = (see Swaters & Balcells 2002). 
^0 takes into account possible deviations from a pure expo- 
nential disc. This may carry valuable information on the 
mass distribution, e.g. if a pseudo-bulge/bar is present, but 
may also refiect variations in the stellar populations and/or 
in the internal extinction, e.g. if the star-formation activ- 
ity is enhanced in the central parts. We use the observed 
central surface brightness ^o. Since we are considering disc- 
dominated galaxies (Sb and later types), we correct /io for 
inclination; we assume an optically-thin disc. Given the am- 
biguity in using either po or /Xd, we include the difference 
A^ = ^d — Mo in the error J^g . This is estimated as 



v/(A/i/2)2 + [2.51og(e)tan(i)5i]2. 



(3) 



For the galaxies from S09, we use the values listed in 
Table A. 5 of Swaters & Balcells (2002) (Harris i?-band). For 
the galaxies from VSOl, we use the surface photometry from 
TuUy et al. (1996) (Cousins i?-band). For the galaxies from 
dB08 and B87, we use the surface photometry from three 
different sources (in order of preference) : Swaters & Balcells 
(2002) (Harris _R-band), Kent (1987) (r-band), and Mufioz- 
Mateos et al. (2009) (Harris _R-band or r'-band). The optical 
filters are comparable, but there can be systematic differ- 
ences of ~0.1 mag (within the typical errors). Two galaxies 
from dB08 (NGC 925 and NGC 7793) and one galaxy from 
B87 (NGC 5371) have no _R-band photometry available and 
have been excluded, reducing the total sample to 52 objects. 

The SO/Sa galaxies from N07 require a different ap- 
proach, because i) the surface brightness rapidly increases 
in the central regions due to the presence of a dominant 
bulge; and ii) several galaxies are at large distances (>30 
Mpc), thus the linear resolution of the optical observations 
is not very high (>150 kpc). For these galaxies, Noordermeer 
& van der Hulst (2007) provide the i?-band disc central sur- 
face brightness /id, extrapolated from an exponential fit and 
corrected for i, and the bulge central surface brightness /ib, 
extrapolated from a Sersic fit to the inner parts after sub- 
tracting the disc contribution. We estimate /io by summing 
the contributions of /id and /ib; the latter value is not cor- 
rected for i as the bulge is assumed to be spherical. The 
errors are given by Eq. 3, where A/i is now the difference 
between /io and the innermost value of /i observed. 



3 THE dRV{Q) -fio SCALING RELATION 

In Fig. 2 (left), we plot no against diiV{0) for the total 
sample of 52 galaxies. There is a clear, striking relation. A 
linear, error-weighted fit to the data yields 



log[dRV{0)] = (-0.22 ± 0.02) no + (6.28 ± 0.40). 



(4) 



As discussed in Sect. 2, the values of dRV{0) for the SO/Sa 
galaxies from N07 are uncertain. However, it is clear that 
these bulge-dominated galaxies follow the same trend de- 
fined by disc-dominated ones. Fig. 2 (left) also shows a linear 
fit excluding the objects from N07 (dashed line). This gives 
only slightly different values of the slope (—0.19 ±0.03) and 
the intersect (5.70±0.57). Considering the different types of 
galaxies and the uncertainties involved, the relation shown 
in Fig. 2 is remarkably tight and extended, spanning more 
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Figure 2. The circular-velocity gradient versus the central surface brightness. The solid and dashed lines show a liner fit to the data- 
points for the total and high-quality samples, respectively. Left: galaxies coded by the reference for the rotation curve (S09; Swaters et al. 
2009, dB08: de Blok et al. 2008, VSOl: Vcrheijen & Sancisi 2001, B87: Begeman 1987, N07: Noordermeer et al. 2007). Right: galaxies 
coded by the value of A/x = fi^ — fiQ (in R mag arcsec"^), that quantifies the "light excess" over an exponential profile. 



than two orders of magnitude in dRV{0) and four orders of 
magnitude in /io. 

Tlie values of the slope and the intersect are likely more 
uncertain than the formal errors, due to several effects in 
the determination of djiV{0) and fiQ. Possible concerns are 
i) the different linear resolutions (in kpc) of the HI and 
optical observations, and ii) the effects on /xq of internal 
extinction, recent star-formation, and/or a LINER core. We 
performed several fits using different methods to estimate /lo 
and duV{0), such as calculating V/R at the innermost point 
of the rotation curve. We obtained slopes always between 
—0.25 and —0.15, and we think that the actual slope must 
be constrained between these values. 

The scatter around the relation is largely due to obser- 
vational uncertainties on diiV{0). Major sources of uncer- 
tainties are i) the galaxy distance, ii) the inclination, and 
iii) the innermost points of the rotation curve (see Eq. 2). 
However, part of the scatter is likely to be intrinsic and 
due to differences in the 3-diniensional (3D) distribution of 
baryons and in the structural component that defines ^o (a 
disc, a bulge, a bar, or a nuclear star cluster). 

To investigate the role played by different structural 
components, in Fig. 2 (right) we plot the same data-points 
coding the galaxies by the value of A/i = /^d — Mo- This 
quantifies the deviation from an exponential law in the inner 
parts of the luminosity profile (in R mag arcsec"^). We dis- 
tinguish between four cases: i) galaxies dominated by a bulge 
(A^ > 1.5), ii) galaxies with a small central concentration 
of light (0.5 $; Afi ^ 1.5) like a pseudo-bulge or a bar, iii) 
galaxies with an exponential disc (—0.5 < A^ < 0.5), and 
iv) galaxies with a central light depression (A^ ^ —0.5). 
The upper-right end of the relation {/iq > 18 mag arcsec"^) 
is populated by bulge-dominated galaxies. It is clear that, 
for these galaxies, the use of ^d instead of /lo would shift 
them away from the relation, as ^d ^ 19-20 R mag arcsec^^ 
(the "Freeman value". Freeman 1970). On the lower-left 



end of the relation, instead, one can find both pure ex- 
ponential discs and galaxies with central light concentra- 
tions/depressions. For these disc-dominated galaxies, the use 
of /id instead of no would still lead to a correlation, but this 
would have a steeper slope (~-0.25). For galaxies with si- 
milar values of ^o, d[iV{0) does not seem to depend on the 
detailed shape of the luminosity profile (simple exponential 
or with a central light depression/concentration). 



4 DISCUSSION 

The correlation between the central surface brightness fiQ 
and the circular-velocity gradient dij1/(0) implies that there 
is a close link between the stellar density and the gravita- 
tional potential in the central parts of galaxies. This holds 
for both HSB and LSB objects, covering a wide range of 
masses and asymptotic velocities (20 < Vnat ^ 300 km s^^). 
The relation between the distribution of light and the 
distribution of mass has been extensively discussed in the 
past (see Sancisi 2004 and references therein). However, 
only few attempts have been made to parametrize this re- 
lation, notably by Swaters et al. (2009) (S09). Fig. 10 of 
S09 plots the logaritmic-slope between 1 and 2 disc scale- 
lenghts Si, 2 = log[V(2/i)/V(/i)]/log(2) versus the "light ex- 
cess" with respect to an exponential disc A/i = /id — Mo ■ It 
shows that a "light excess" (a bulge-like component) cor- 
responds to a "velocity excess" in the rotation curve with 
respect to the expectations for the underlying exponential 
disc. The limitation of that parametrization is that it does 
not capture the dynamical difference between HSB and LSB 
discs, that are known to have steeply-rising and slowly-rising 
rotation curves, respectively (e.g. TuUy & Verheijen 1997). 
In Fig. 10 of S09, indeed, both HSB and LSB exponential 
discs have A/i ~ and Si, 2 — 0.5. The latter result is due to 
the fact that Si, 2 is, by definition, a scale-invariant quantity. 



© 0000 RAS, MNRAS 000, 000-000 



A scaling-relation for disc galaxies 5 



that is not expected to depend on /^o or Knax- In contrast, 
rfH^(O) measures the inner slope of the rotation curve in 
physical units (km s~^ kpc~^) and is directly related to the 
central dynamical surface density (in Mq pc~^), providing 
insights in the underlying physics, as we now discuss. 

For a 3D distribution of mass, the rotation velocity V of 
a test particle at radius R is given, to a first approximation, 
by 



R 



GM, 



dyn 



i?2 



where G is Newton's constant, M^ 



(5) 
4/37r7?3p is the 



Jdyn */ >^ rdyn 

dynamical mass within R, and a is a factor that depends on 
the detailed mass distribution (for a spherical distribution 
of mass Q=l, while for a thin exponential disk a ~ 0.76 at 
R = 0.5Ro). For i? -> 0, we have 



dV V 



^ = ;^ = y^G^yn.o=^/3G^^^^^^^ 



,/Obar,0 



(6) 



where /3 = A/Sna, pdyn,o and pbar.o are, respectively, the 
central dynamical and baryonic mass densities, and /bar.o = 
Pbar,o/Pdyn,o IS the baryou fraction in the central regions. 
Note that /bar.o niay strongly differ from the "cosmic" 
baryon fraction, and can vary widely from galaxy to galaxy, 
depending on the formation and evolution history. Observa- 
tionally, we measure po which is related to pbar.o by 



^0 = -2.51og[pbar,() A2:(Mbar/i)' 



(7) 



where Az is the typical thickness of the stellar component 
(either a disk or a bulge) and Mbar/i is the baryonic mass- 
to-light ratio, including molecules and other dark baryonic 
components. Thus, we expect the following relation 



log[dfly(0)] = -0.2 po + 0.5 log /3G 



Mbar/L 



Az/ba 



(8) 



In Sect. 3, we mentioned that the slope of our rela- 
tion is not well-determined due to several uncertainties in 
the measurements of dRV{0) and po- However, it is con- 
sistent with —0.2 and can be constrained between —0.15 
and —0.25. Were the slope exactly —0.2, the second term of 
Eq. 8 would be a constant, implying a puzzling fine-tuning 
between the 3D distribution of baryons (/3 and Az), the 
baryonic mass-to-light ratio (Afbar/i), and the dark matter 
content (/bar.o)- 

Despite the uncertain value of the slope, the results pre- 
sented here show a clear relation between the central stellar 
density in a galaxy and the steepness of the potential well 
(see also Sancisi 2004; Swaters et al. 2011). This implies 
a close link between the density of the baryons, regulated 
by gas accretion, star-formation, and feedback mechanisms, 
and the central density of the DM halo, together shaping 
the inner potential well. This may represent a challenge for 
models of galaxy formation and evolution. Future observa- 
tional studies may help to better constrain the slope of the 
relation, while theoretical work should aim to understand 
its origin. 



5 CONCLUSIONS 

We measured the circular-velocity gradient diiV{0) for a 
sample of spiral and irregular galaxies with high-quality 



rotation curves. We found a linear relation between 
log[dal/(0)] and the central surface brightness po with 
a slope of about —0.2. This is a scaling-relation for disc 
galaxies that holds for objects of very different morpholo- 
gies, luminosities, and sizes, ranging from dwarf irregulars 
to bulge-dominated spirals. This relation quantifies the 
coupling between visible and dynamical mass in the central 
parts of galaxies, and shows that the central stellar density 
closely relates to the inner shape of the potential well. 
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